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AN INVESTIGATION OF SELF-SUBTRACTION 
HOLOGRAPHY IN LiNbO^ 

D. W. Vahey, R. P. Kenan, N. F. Hartman, anf5 R. C. .Sherman 
Battelle Columbus Laboratories 

SUMMARY 

This report describes the results of a fifteen-month investigation 
of self-subtraction holography in LiNbO^ optical waveguides. Self subtraction 
is the term used to describe the gradual transfer of energy from one optical 
beam to another that can occur when the two beams are used to record a holo- 
gram in LiNbO^ or other photorefractive materials. The effect is potentially 
useful in optical processors where the information signal is superimposed on a 
slowly varying and uninteresting background signal. Self subtraction affords 
the possibility of eliminating the background from the total signal received. 

The quality of self subtraction is inferred from measurements of the 
time required to form the hologram, the minimum beam intensity generated, the 
ability of the LiNbO^ to form many successive holograms without fatigue, and 
the stability of the holograms during long-term readout. In this program we 
tested in depth one sample having self-subtraction characteristics that were 
very promising: hologram-formation times were on the order of 150 sec, the null 
signal was less than 2.5% of the peak signal, and no fatigue nor instability 
was detected over the span of the experiments. Another sample, fabricated 
similarly or with slight modifications did not perform nearly as well. In all 
samples, attempts to improve self-subtraction characteristics by various 
thermal treatments. had no effect or had adverse effects, with one exception in 
which improvement was noted after a time delay of several days. 

A theory developed to describe self subtraction showed the observed 
decrease in beam intensity with time, but the shape of the predicted decay curve 
was oscillatory in contrast to the exponential-like decay observed. The theory 
was also inadequate to account for the experimental sensitivity of self sub- 
traction to the Bragg angle of the hologram. 



We conclude that self subtraction is a viable method for optical pro- 
cessing systems requiring background discrimination. However, understanding 
the physical basis for the effect has only begun and further work of both an 
experimental and theoretical nature is required in order to optimize and con- 
trol self subtraction for use in applications. 
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I . INTRODUCTION 


DEFINITION OF SELF-SUBTRACTION HOLOGRAPHY 


Self-subtraction holography is the process by which one of two beams 
used to write a hologram gradually transfers its energy to the other beam 
during hologram formation, and emerges from the hologram region with decreasing 
intensity. The designation of the term "self-subtraction" to describe this 
phenomenon is not obvious unless one views it as a special case of the tech- 
nique known as holographic subtraction. In this technique, shown in Fig. 1, 
a viewer or other sensing system detects the coherent superposition of two 
wavefronts. One is the instantaneous wavefront from a specific object that 
passes through a hologram of that object without being diffracted. The 
other wavefront is generated by reconstruction in the hologram. If the 
relative phase of the two wavefronts is adjusted to tt, they will destructively 
interfere and the viewer will see a dark field. It is said that the re- 
constructed wavefront is subtracted from the real wavefront; hence we have 
the terminology "holographic subtraction"^^^ In this example it was implied to 
be necessary to adjust the relative phase shift between wavefronts to tt by 
an external means, such as a phase shifter. In certain holographic materials, 
such as LiNbO^, the relative phase between the reconstructed wavefront and 

the real wavefront is equal or close to tt by virtue of the hologram-formation 

( 2 ) 

mechanism. In this case, an externally applied phase adjustment is not 
necessary to achieve the subtracted condition. The situation here is there- 
fore called self-subtraction. 

ADVANTAGES OF SELF-SUBTRACTION HOLOGRAPHY 

The potential applications for self-subtraction holography are the 
same as those for subtraction holography in general. By use of this tech- 
nique, one hopes to more readily detect small changes in objects or scenes 
under study. A specific application which has been a motivating force for 
this investigation is that of the Integrated Optical Data Preprocessor, 
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Phose Shifter 



Hologram being 
written 


Phase object 


(a) THE FORMATION OF THE HOLOGRAM OF A PHASE OBJECT 




PHASE SHIFT WILL RESULT IN THE OBSERVER SEEING 
A DARK FIELD 


FIGURE I-l. the HOLOGRAPHIC SUBTRACTION PROCESS. 
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currently under development at Battelle.'- In this application the "object" 
is a spatial phase modulation imposed on an optical guided wavefront. The 
phase modulation is achieved by applying voltages from an array of sensors to' 
an array of electrodes spanning a waveguided beam in electrooptic LiNbO^- 
Changes in the wavefront are detected by holographic subtraction and provide 
a direct indication of changes in the voltages provided by the sensors. The 
usual situation is that in which sudden changes in sensor voltages represent 
phenomena that are of interest, while slow variations reflect changes in the 
ambient background and are to be ignored. Conventional holographic subract- 
tion, in addition to requiring an externally applied phase shift, uses per- 
manent holograms geared to specific sensor outputs, and so indiscriminantely 
detects both fast and slow changes in these outputs. However, self-subtraction 
holograms are continuously updated in time, and so detect only those sensor out- 
puts that change more rapidly than the characteristic update time constant, as 
required. 


ACCOMPLISHMENTS OF THE SELF-SUBTRACTION PROGRAM 

With this system application in mind, it is appropriate to describe 
the accomplishments of the self-subtraction program. Working with holograms 
formed in LiNbO^ optical waveguides, we performed experiments which served to 
characterize the phenomenon with regard to (1) the extent to which the sub- 
tracted optical beam could be extinguished, (2) the time constants associated 
with hologram formation, and (3) the material fatigue associated with the con- 
tinuous updating of the self-subtraction hologram. Our results in the charac- 
terization phase were very encouraging, showing 40-to-l maximum-to-minimum 

2 

signal ratios, time constants on the order of 10 sec, hologram stability great- 
er than 7 hours under contant exposure, and hologram recycleability greater 
than 9 cycles without any signs of degradation. 

We also performed experiments designed to elucidate the physical 
mechanisms for self-subtraction, and to optimize the effect for use in holo- 
graphic optical-processing systems. These experiments involved the modification 
of optical beam ratios, beam intensities and beam widths, the use of both Ti- 
diffused and Li-outdif fused LiNbO^ waveguides, the variation of the hologram 
grating period, and the use of various reducing/oxidizing thermal treatments. 
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Although a variety of interesting phenomena were observed and a reasonably re- 
producible procedure was determined for obtaining useful self-subtraction 
characteristics in LiNbO^ waveguides, these experiments were not particularly' 
useful in elucidating the physical mechanisms of self subtraction. A theory 
was developed based on the simplest physical model that appeared to have the 
sophistication necessary to describe the basic self-subtraction phenomenon. In 
so far as it was possible to test the theory during the program, however, the 
agreement with experiment left something to be desired. 

In summary, we were successful in characterizing the phenomena of 
self subtraction experimentally, and in determining a set of conditions for 
producing good self subtraction in outdif fused LiNbO^ waveguides. The results 
were encouraging with regard to the eventual use of the phenomenon in practical 
systems. Yet much work remains to be done to understand the physics of self 
subtraction and to optimize the effect for use in applications. 
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II . TECHNICAL BACKGROUND 


In this section we describe our current understanding of the physical 
basis for self-subtraction holography. First we describe the spatial phase re- 
lationship between the index-modulation pattern of a hologram and the interfer- 
ence pattern of the writing beams that leads to subtraction. Then we describe 
how this spatial relationship may result automatically in the case of photo- 
refractive holographic materials like LiNbO^. Then we develop a more rigorous 
description of self subtraction based on a common model for photorefractivity, 
and we indicate some of the behavior predicted by this model. This will lead 
us naturally into a discussion of our experiments and experimental results in 
Sec. III. 


A PHYSICAL DESCRIPTION OF SELF "SUBTRACTION 


Figure II-l shows the interference pattern formed by two beams of 
unit intensity that have combined to form a hologram in a photosensitive mate- 
rial, The drawing indicates schematically that the relative phase between the 
interference maxima and the refractive index maxima of the hologram is zero; 
that is, intensity maxima and index maxima overlap. The symmetry of the situa- 
tion dictates that the output intensities of the two optical beams will each be 
unity, since they were taken to be unity at the hologram input. In any event, 
the output beams may be regarded as being composed of diffracted and undif- 
fracted components of the input beams R^ and S^. If n is the diffraction effi- 
ciency of the hologram, we may write 


S . 

R - + e’-'*’ 




(II-l) 


where the first term on the right of each equation is the undiffracted contri- 
bution to the output amplitude, and the second is the diffracted contribution. 
A natural phase shift has been associated with the diffraction process. The 
existence of such a phase shift is seen to be imposed by the symmetry 
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FIGURE II-l. INTERFERENCE PATTERN AND CORRESPONDING REFRACTIVE-INDEX 
PATTERN AFTER HOLOGRAM FORMATION. INDEX MAXIMA ARE 
ALIGNED WITH INTENSITY MAXIMA (NO PHASE SHIFT) . 



2 2 2 2 
requirement that S = R =1. (Recall we assumed = R^ = 1 at the hologram 

input.) This condition is satisfied by Eq. (1) only if = + ir/2. The cor- 

rect choice 


ip = -ir/2 

(II-2) 

(4) 

has been derived by Staebler and Amodei. 


Let us now apply an external phase shift if> to 

the signal beam in 

Fig. II-l, so that S^ = 1 e^*^ while R^ = 1 as before. 

Equation (1) becomes 

s - (i-n)"/' 

(n-3) 




The output intensities are 

SS* = 1 - 2(n)^^^ (1-n)^^^ sin (j) (II-4) 

RR* = 1 + sin <l> . 

When <(> = ir/2, the signal beam reaches its minimum value and may be described 
as being holographically "subtracted". In contrast, the reference beam is 
said to be holographically "added". 

Figure II-2 shows the new spatial relationship between the hologram 
index pattern and the interference pattern. It suggests that in the absence of 
an external phase shift (j) = ir/2, subtraction of the signal beam could have been 
obtained by displacing the hologram pattern by 1/4 fringe, as this accomplishes 
the same quadrature phase relationship shown in the figure. Materials like 
LiNbO^ for which this phase relationship can be established naturally through 
the hologram-forming mechanism are those which exhibit self subtraction. In 
LiNbO^, the dominant hologram-formation mechanism is termed photorefraction, 
which we now discuss. 


THE PHOTOREFRACTIVE PROCESS 


The term photorefraction designates photoinduced refractive- index 
changes. In LiNbO^ the refractive- index changes arise from the generation of 
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FIGURE II-2. INTERFERENCE PATTERN AND INDEX PATTERN AFTER 

A tt/2 phase SHIFT APPLIED TO THE SIGNAL BEAM. 


a nonuniform distribution of electrons trapped at impurity sites within the 
crystal. When the crystal is illuminated, electrons are excited into the 
conduction band where they drift and diffuse until they are retrapped. For 
a spatially structured light beam, such as an interference pattern, the general 
tendency is for the population of trapped electrons to increase in the dark 
regions of the crystal and to decrease in the illximinated regions. The process 
is shown schematically in Fig. II-3. 

The associated space-charge electric fields produce a refractive- 
index change by means of the linear electrooptic effect. The process continues 
uniformly until the space-charge fields become so large that they oppose fur- 
ther migration of electrons . Often the fields are as large as 10^ V/cm, 

-4 

corresponding to refractive- index changes greater than 10 . These index 

changes are large enough to produce holograms having diffraction efficiencies 
approaching unity. 


The Origin of Donors and Traps 

Photorefractive sensitivity is observed to be largest in LiNbO- 

(5) 

crystals that have a high concentration of iron. Since iron is a naturally 

occurring impurity in LiNbO-, even samples that are described as nominally 

^ 2 + 
pure have sufficient ifon for the formation of efficient holograms. Fe im- 

3_l_ ( 8 ) 

purities are known electron donors, and Fe"^ impurities are known traps. 

However, photorefractive effects have also been observed in crystals doped 

(9) 

with other ions, such as Cu and Mn, and holograms formed by two-photon 
processes appear to rely on transient photoinduced color centers for sensitiza- 
tion. Trapping sites necessary for the formation of long-duration holo- 

grams may be associated with crystal imperfections other than impurities, such 
as oxygen-site vacancies. 

The Migration of Photoexcited Electrons 

After an electron is excited into the conduction band, it may be 
retrapped at the same location or it may migrate to other regions of the 
crystal. This migration is partially diffusion controlled, but a more signifi- 
cant mechanism is electron drift associated with a bulk photovoltaic 
effect. This effect originates in the pyroelectric nature of the LiNbO^ 
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FIGURE II-3. CONCEPTUAL DRAWING SHOWING REDISTRIBUTION OF 
FILLED TRAPS DURING ILLUMINATION BY AN INTER- 
FERENCE PATTERN. OPEN CIRCLES ARE UNFILLED 
TRAPS, CLOSED CIRCLES ARE FILLED TRAPS. SOLID 
CURVE DESIGNATES INTENSITY VARIATION. DASHED 
CURVE DESIGNATES REFRACTIVE INDEX VARIATION. 





host. Fe impurities are separated from neighboring Nb ions by distances that 
vary in the positive and negative c-directions. Electrons excited from Fe 
impurities consequently exhibit a preference to motion in one of the two 
directions, along which they drift uniformly until they are retrapped. The 
migration establishes a photovoltage across the illuminated region, the sign 
of which shows that electron migration occurs along the - c axis. The effect 
is similar to that which would be observed if the electrons drifted in response 
to a positively directed electric field. Indeed, the migration can be speeded 
or slowed by applying a suitable voltage to the crystal. The field required 
to counteract the photovoltaic effect is found to be about 5 x 10^ V/cm.^^^^ 


The Photorefractive Equations 

At any point within the crystal, the concentration of conduction 
electrons n can be described by a continuity equation. Optical excitation 
generates conduction electrons, while retrapping, photovoltaic drift, and 
diffusion remove them. Kim et al^^^ have employed an equation which appears 
to account for all relevant processes: 

[3/3t+l/T-D3^/3x^+y3/3x(E+E(x))]n = I(x)Naq/hv (II-5) 


The terms in this equation have the following meaning: 

X = the direction parallel to the optic axis; 

3n/3t = rate of change of electron concentration; 

n/x = rate of electron retrapping; 

2 2 

D3 n/3x = spatial derivative of the diffusion current, D3n/3x, where 
D is the diffusion coefficient; 

y3/3x(En+En)= spatial derivative of the drift current, where E is the effec- 
tive photovoltaic field, E is the space charge field, and y is 
the electron mobility; 

I(x)Noq/hv = rate of generation of conduction electrons per unit volume, 

where N is the trapped electron concentration, a is the absorp- 
tion cross section, q is the quantum efficiency for processes 
resulting in migration, hv = photon energy, and I(x) = intensity. 
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For simplicity we eliminate diffusion from further consideration be- 
cause holographic sensitivity in LiNbO^ is much greater for geometries in 
which both drift and diffusion operate than for geometries in which only dif- 
fusion operates. We also neglect the effects of E on electron drift, since 
it is small during the initial stages of hologram formation in which we have 
interest. Finally, we neglect 8/9t in comparison to 1/x, since the retrapping 
time is expected to be on the order of picoseconds while changes in n from the 
holographic process occur over a span of seconds. This simplifies Eq. (5) to 

n + L(9n/3x) = I(x)Naqx/hv 5 g(x)x (II-6) 

where L = pEx is the characteristic length that a conduction electron drifts 
before it is retrapped, and g(x) is a simplified notation for the generation 
rate. In the case of illumination by the interference pattern 

g(x) = §q( 1 + ® 2kxsin0) , (II-7) 


the solution to Eq. (6) is 

n(x) = n^ + n^ sin(2ksin6) + n^ cos(2kxsin0) , (II-8) 

So^ ’ 

g^x m X/(l + X^) 
g^x m/(l + X^) 

2kLsin0 , 

One particular interest is in the concentration of trapped space charge 
p described by the continuity equation p = eyE9n/9x. After determining p with 
the help of Eq. (8), Maxwell's equation 9E/9x = p/c is integrated to obtain 
the space-charge electric field. The associated refractive index perturbation 
is n^ = (1/2) n^r^^Ej where n^ = 2.2 is the extraordinary refractive index of 
LiNbO^j and r^^ = 30 x 10 m/V is the appropriate electrooptic coefficient. 

We find 
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(II-9) 


cos(2kxsin0 - (j)^) 

3 

n^g = t(n^rgg/2) (eLg^/e)m coscj)^ 

(i) = tan ^ 2KLsin0 . 

g 

The phase shift between the hologram index pattern and the optical 

interference pattern is seen to arise as a natural consequence of electron 

drift. The value (j> = tt/2 required for self-subtraction can be obtained by a 

g 

crystal for which 2kLsin0>>l; that is, the drift length is long compared to 
the hologram fringe spacing. 

The value of L can be experimentally increased by reducing the num- 
ber of empty traps. This can be done by increasing the relative number of 
filled traps, using reducing treatments of the type described by Phillips and 
Staebler.'' ^ This was attempted during the program and the results are re- 
ported in the next section. 

It is worth noting that in crystals for which electron diffusion is 

the dominant transport mechanism a phase shift <t> = tt/2 automatically occurs. 

g 

This may be seen by keeping D and eliminating y in Eq.' (5), and following pro- 
cedures similar to those we have already used in analyzing that equation. 

THEORETICAL ANALYSIS OF SELF -SUBTRACTION 


In this subsection we present a theoretical analysis of grating 
formation in photorefractive LiNbO^ for the interesting case of a dominant 
migration path length L. The treatment is a coupled-wave analysis of Maxwell's 
equations as they are constrained by the photorefractive equations of the pre- 
vious subsection. Time-dependent Maxwell's equations are used in view of the 
slowly varying temporal changes of holograms during their recording. A treat- 
ment intermediate to the qualitative picture of self subtraction that intro- 
duced the section and the more rigorous derivation presented here is contained 
in Appendix A. 

We assume two coherent beams propagating at angles + 0„ relative to 
the Z-axis of the LiNb 03 material, as shown in Fig, II-l, These are the 


15 



reference and signal beams, R and S described by amplitudes R and S, phases 
(j)j^ and (J)g, and propagation constant g as follows: 


R = R(z,t) giScosez + igsinex 


(II-IO) 


S = S(z t) gi3cos0z - iBsin0x 


Note that previously R and S were complex and played the role now played by 
R and S. From this point forward, R and S are real quantities. For simplicity 
the beams are assumed to be polarized in the y direction. The total electric 
field E = R + S satisfies the wave equation 


V^E + k^E = 0 (II-ll) 

where k = 2im(x, z, t) /X^, n(x,z,t) is the refractive index and is the free- 
space wavelength. The spatial and temporal variation of the refractive index 
reflects the fact that the field is writing a hologram by means of the photo- 
refractive characteristics of the material. Generalizing the result of 
Eq. (9), we assume a refractive index variation of the form 


n(x,z,t) = + n^Q(z,t)cos[2kxsin0g+(|>j^(z,t)-({ig(z,t)-(i)g(z,t) ]. (11-12) 


Here n^ is the average refractive index. Accordingly, the propagation constant 
6 of Eq. II-IO is given by S = 2-!rn^/X^. <j>g(z,t) describes the Instantaneous 

phase shift between the hologram index pattern and the optical Interference 
pattern. 

2 2 2 

By substituting k = k^n (z,t) into the wave equation along with 
E = R + 5, we can reduce it to coupled equations for the individual beam am- 
plitudes R and S. The procedure for this is exactly that originally employed 
( 12 ) 

by Kogelnik. Assuming n^^Q<<n^ and neglecting all but synchronous terms 

in the wave equation, we find 


_• (<l> +<{>c5-0 

3R/3z - i(k n,„/2cos0„) e ° S = 0 

O XU D 


i(4>g+<)>s-<^R) 


R - 0 . 


(11-13) 
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3S/3z - i(k nT«/2cos0„) e 

O iU D 



To complete the formal description of grating formation, it is neces- 
sary to develop equations that display the dependence of the grating parameters 
n^Q and <{>g the fields R and S. The optical interference pattern at any 
value of z and t is 

I = n^(e^/4y^)^'^^{R^(z,t) + S^(z,t) + 2R(z,t)S(z,t) (11-14) 

cos[2gxsin9g - <j)g(z,t) + ij)^(z,t)]} 

This interference pattern excites trapped electrons into the LiNbO^ conduction 
band where they migrate antiparallel to the crystal's optical axis as described 
above. Consider electrons initially in a volume element AV = £, AxAz, where I 

y y 

is the length of the hologram in the y direction. The number of electrons ex- 
cited in an interval At is called AN_, where 

AN_(x,z,t) = (N^JlyAxAz)F^aAtn [1 + m cos(2gxsin6g-(t)g+(|)j^) ] . (11-15) 

In this equation is the density of trapped electrons, 

Fo = [n^(c^/4y)^''^ (R^ + S^)]/hv (11-16) 

is the average flux density of the beam, hv being the energy per photon, a is 
the absorption cross section, n is the quantum efficiency for the generation of 
conduction electrons, and 


m = n (c /4y)^'^^ 2RS/F hv (11-17) 

o o o 

is the visibility of the interference pattern. The result of photorefractive 
migration is to distribute these electrons into new volume elements H^AxAz 
centered at x' > x. The fact that these electrons are lost to the original 
volume element is indicated by use of the notation AN_. However, these elec- 
trons are replaced by electrons generated at x' < x which drift into the volume 
element at x. This contribution is 

AN^(x,z,t) = I ^ AN_(x',z,t) (11-18) 
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where L is the migration path length, and the exponential term accounts for 
the decreasing contribution from increasingly distant volume elements. 

The total change in the number of trapped electrons at x is 
AN = AN_|_ - AN . An expression for AN is obtained by carrying out the inte- 
gral of Eq. (18) and by subtracting from the result the expression of Eq.(15). 
We find 

AN= N (5, AxAz)F apAt m(l + KL (II-19) 

o y o 

{-KLcos(Kx + ([))+ sin(Kx + i^)) } , 

where we have introduced the grating vector magnitude K E 26sin9„ and the 
phase difference (fi E - (j)^ . 

The redistribution of electrons described by Eq. (19) produces an 
incremental space-charge electric field AE consistent with the Maxwell equa- 
tion 

3(AE)/3x = Ap/e^K , (11-20) 

Ap = -eAN/ (JlyAxAz) , 

where Ap is the space charge density and k is the dielectric constant. The 
space-charge field AE in turn causes a proportional refractive index change 
of electrooptic origin, An, given by 

An = ^ n^r AE , (11-21) 

z o 

where r is the appropriate electrooptic coefficient. By inserting Eq. (19) 
into Eq. (20), integrating, and then substituting the result for AE into 
Eq. (21) we find 

-j P P -1 /2 

An(x,z,t) = At(n\/2)(eL/e^K)N^F^man(l + kV) ‘ (11-22) 

cos(Kx + <f) - Q) , 
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where Q = tan ^ KL is the photoref ractive phase shift obtained previously. 

The term F^m on the right-hand side of this equation is related to the field 
amplitudes R and S by Eq. (17), which we restate: 

F^m = (hv) ^ R(z,t)S(z,t) (11-23) 

A differential equation relating the refractive index of the material to the 
fields is obtained by inserting this result into Eq. (22), dividing by At, and 
letting At become infinitesimal. We find 

3n'(x, z, t) /3t = A R(z,t)S(z,t)cos(Kx + -Q) (11-24) 

K b 

A = (n^r/2) (eL/ic) (X /h)N ori(l + . 

o o o 

It is convenient to rewrite this result in the form 

3n/3t = ARS cos(Q-(J>)cosKx + ARS sin(Q-cj)) sinKx . (11-25) 

For comparison, Eq. (12) defining n(x,z,t) is rewritten 

n = n^ + n^Q cos(t()g-^)cosKx + n^Qsin({j>g-(j)) sinKx . (11-26) 

By taking the time derivative of this expression and comparing the coefficients 
of the cosKx and sin Kx terms that result with those of Eq. (25), we can make 
the identification 


3/3t (n^Qe^^‘*’g ‘*’^) = ARS e^^'^ . (II-27) 

This equation together with Eq. (13) provides the information we require to 
obtain the field and refractive- index variations encountered during photo- 
refractive grating formation. Analysis of the three complex equations is 
facilitated by rewriting them as six real equations for the quantities R,S, 
(})^, (J>g, n^Q and (J)^. This is straightforward to carry out. The resulting 
equations are 
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R 

z 

- BSn,rtSiniJ) 
10 ^g 

= 0 

(a) 

(11-28) 

S 

z 

+ BRn-»sindi 

10 ^g 

= 0 

(b) 



- BSn, _cos(j> 

10 ^g 

= 0 

(c) 


- 

- BRn, -cosij) 

10 ^g 

= 0 

(d) 



(n^Q)j. - ARScos(Q-<J)g) = 0 (e) 

nio(<()g-<(>)t - ARSsln(Q-<j)g) = 0 (f) 

We have introduced B = k /2cos9„, and we have employed subscripts z and t to 

O o 

designate 3/3z and 3/3t, respectively. The initial conditions and boundary 
conditions to use in solving these equations are as follows: 

R(0,t) = R(z,0) = R^ (11-29) 

S(0,t) = S(z,0) = 

(})j^(0,t) = (|>j^(z,0) = 0 

(j)g(0,t) = lj)g(Z,0) = 0 

<{>g(0,t) = <))g(z,0) = Q 

nj^g(z,0) = 0 

In general, these equations must be solved numerically. There are, 
however, two special cases for which analytic solutions are available. These 
are the cases Q = 0 and ir/2. In the former, we also require R^ = S^. In 
the latter, we require R^>>S^. We now consider these special cases, and, in 
closing, the general case. 

Case I. Q = 0, R^ = S^ . With Q = (j)g(z,o) = 0, Eqs. (28a) and (28b) 
show that there is no tendency for R and S to deviate from their initial values 
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R and S , at least for short writing times. The result is that cj) and ({) in- 
0 0 o 

crease in a parallel manner [Eqs. (28c) and (23d) ] with the difference (j) = 0 

maintained. This forces 6 to remain at zero in accord with Eq. (28f) . The- 

§ 

same argument can be employed to describe successive small increments of 
writing time. The net result is 


R(z,t) = R^ (11-30) 

S(z,t) = 

<j)j^(Z,t)-$g(Z,t) = 0 

nio(Z’t) = AR^S;t 
<|)g(z,t) = 0 . 


As time increases, the space-charge fields eventually assume a major role in 
governing the redistribution of photoexcited carriers. This has not been taken 
into account, and is the source of the unphysical result that n^^Q is propor- 
tional to t at all writing times. However, high diffraction efficiency grat- 
ings can be recorded in times for which Eq. (30) is valid. 

Case II Q = ir/l, R^»S^ . With Q = (j) (z,0) = ir/2, Eqs. (28c) and 

(28d) indicate no initial change in c|) or (}>£,. This is true even though there 

K b 

is a change in R and S, with R increasing from R^, and S decreasing from S^. 

By Eq. (28f) , the stability of <|)„-({)„ insures the continuation of the initial 

K b 

condition <p = tt/ 2, so that the same argument continues to apply for all writ- 
8 

ing times. The coupled Eqs. (28) simplify to 


(n 


BSnio - 

0 

(a') 

(H-31) 

+ BRn^^Q = 

0 

(b') 


- ars = 

0 

(e*) 



To obtain a single partial differential equation for S, we make use of the 
energy-conservation condition 
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(11-32) 


2 2 2 2 2 
R + S = R + S = U 
o o o 

which may be derived from Eqs. (31a*) and (30b*). By substituting 
R = (U^ - in Eq. (31b*), dividing by (U^ - to segregate terms 

containing S from terms containing differentiating with respect to time, 

and rearranging terms, we obtain 

(U^-S^)S , + SS + AB S(U^-S^)^ = 0 . (11-33) 

ozt zt o 

This may be converted to an ordinary differential equation by introducing the 
variable w = ABzt. The result is 

(wS + S )(U^ - S^) + wSS^ + S(U^ - S^)^ = 0 . (11-34) 

ww wo wo 

This equation is valid for all initial values of U and S . However an analytic 

2 2 o o 

solution is obtained when U^>>S , that is, when most of the energy is confined 
to the reference beam. On the basis of this assumption we neglect terms that are 
quadratic in S or higher in Eq. (34), with the resulting simplification: 

wS + S + S = 0 . (11-35) 

ww w 

(13) 

This is a member of the class of equations solved by Bessel functions.'' ^ The 
solution is 

S(z,t) = S^J^[2(ABzt)^^^] . (11-36) 

2 

An example of the shape of the curve (S/S^) is indicated in Figs. II-4 and 
II-5 along with numerically generated curves obtained from the computer solu- 
tion of the general equations, Eq. (28) » which we now discuss. 

Case III. The General Case . A FORTRAN computer program was develop- 
ed to solve Eq. (28) . The numerical method employed was that of finite differ- 
ences, modified to insure conservation of energy of reference and signal beams. 
The program has not been optimized. Rather precision of results has been 
achieved by the procedure of using sufficiently small increments of space and 
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NORMALIZED SIGNAL BEAM INTENSITY, S (1 cm,t)/S 



FIGURE II-4. SIGNAL-BEAM INTENSITY VS RECORDING TIME. SOLID CURVE: 
COMPUTED FOR (Sq-Rq)^ = 0.01, Q = 1.5. CROSSES: AP- 

PROXIMATE BESSEL-FUNCTION SOLUTION. DASHED CURVE: COM- 
PUTED WITH GREATER PRECISION. 
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NORMALIZED SIGNAL-BEAM INTENSITY, S (1 cm,t)/(R' 
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FIGURE II-5. SIGNAL-BEAM INTENSITY VS RECORDING TIME. SOLID CURVE: 
COMPUTED FOR s2 = r2 = 1/2, Q = 1.4. DASHED CURVE; 
BESSEL FUNCTION APPROXIMATION. 



time to integrate the coupled equations. The difficulty in implementing a pro- 
gram to solve the given system may be inferred from an examination of Eq. (28d) . 
As the S-beam amplitude decreases toward zero during subtraction, the spatial- 
rate of change of the phase becomes very large. This jeopardizes the preci- 
sion of the calculation unless care is taken. Problems that we encountered in 
this regard prevented us from performing a very detailed study of the implica- 
tions of Eq. (28). We have, however, identified two important features: 

(1) Slight reduction of the parameter Q from tt/ 2 does not seriously degrade 
the crystal's capability for performing self subtraction; and (2), self sub- 
traction is largely independent of the initial beam ratio (S^/R^) . 

Figure II-4 shows the computer generated self-subtraction curve for 
the case (S^/R^) = 10 and Q = 1.5. The value for A was 10 sec ; and for 

B, 4.971 ym Crosses placed close to the computed curve show the correspond- 
ing value of the Bessel-function approximation given by Eq. (36). There is 
some evidence that the agreement between the two calculations is even better 
than shown. Over a limited range that was tested, the computed curve approach- 
es the Bessel-function curve as more precision is obtained from the calcula- 
tion (dashed curve in Fig. II-4) . 

In Fig. II-5, we present the computer-generated results for the case 

2 

(S^/R^) = 1 and Q = 1.4. We find, first of all, a secondary maximtim that is 

about 20% of the initial intensity. This was also the case for the beam ratio 

2 -2 

(S /R ) =10 . The nature of the subtraction curve does not, therefore 

appear to depend strongly on beam ratio. One noticeable difference between 
the computed curves of Figs. II-4 and 5, however, is the extent to which 
minima approach the baseline. The approach is closer in the case of the smaller 
beam ratio of Fig. II-4, which incidentally pertains to a value of Q that is 
closer to the ideal subtraction value x/2. For the curve of Fig. II-5, the 
second minimum is significantly further from the baseline than the first minimum. 

This could be an indication that the long-term signal intensity does not ap- 
proach zero as the Bessel-function approximation suggests. We looked for this 
Instability experimentally but did not find it. Details of the experiment will 
be provided in Sec. IV. It is possible that photorefractive saturation effects 
not accounted for in the present theory could prevent instabilities in the sub- 
traction process from becoming evident, thus explaining our experimental result. 
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Figure II-5 shows, in addition to the computed curve, the Bessel-func- 
tion substraction curve obtained from the theory for the case of a small beam 
ratio and Q = ir/2. At small times, differences in the computed and Bessel 
curves reflect differences in the beam ratio. As subtraction progresses, the 
beam ratio becomes small even though it was not small initially. We therefore 
expect and observe oscillatory behavior similar to that shown by the Bessel- 
function curve of Fig. II-5. Failure of the computed curve to reach the base- 
line is probably more an indication of the choice Q = 1.4, that is, Q less than 
tt/ 2, than it is an indication of beam ratio. If this is true and if the failure 
of the computed curve to reach the baseline is an indication of self-subtrac- 
tion instability, then the time required for instability to become evident 
probably decreases as the deviation of Q from the ideal subtraction condition, 

Q = it/2, increases. This is a question to be addressed in future theoretical 
work on self subtraction, along with changes in the theoretical model to account 
for photorefractive saturation effects that could limit the appearance of in- 
stabilities and smooth the oscillations of the subtraction curves. 

Thus far our treatment of the general Eq. (28) has addressed only the 
evolution of beam intensity during the subtraction process. In concluding this 

section we discuss the beam phase 4) = ^ the limiting case describ- 

Kb b 

ed by the solution 

1 /2 

S(z,t)/S^ = J^[2(ABzt) ' ] , 11-37 

zero crossings of the Bessel function correspond to points where the phase of 
the signal beam, (jigi jumps by r. In the general case the change in phase is 
continuous, but rapid variations occur at times when the signal beam intensity 
is at a minimum. Figure II-6 and II-7 show the computed variations of the 

phase 4>c(l cm,t) for the initial conditions used to generate Figs. II-4 and 

^ 2 
II-5, respectively. Arrows designate those times at which S is a minimum. 

In the case of Fig. II-6, for S^/R^ = 1/10 and Q = 1.5, the off-scale values 

of <)>g near the first minimum of result from imprecise computation in that 

region. A dashed curve shows more precise values. 
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SIGNAL BEAM PHASE, - 4 . (1 cm,t) 



SIGNAL-BEAM PHASE, -4>„(1 cm,t) 





STATUS OF THE THEORETICAL UNDERSTANDING OF 
HOLOGRAM FORMATION IN LlNbO ^ 

The level of the theory presented in this section is the simplest 
that one can appeal to and still describe beam- interaction effects such as 
those associated with self subtraction. This is evidenced by the fact that the 
very complicated photorefractive effect is reduced in this work to a single 
parameter L, the electron drift path length. This was appropriate because any 
attempt to fit experimental results to a theoretical model should begin with 
the simplest model possible. The reader should understand, however, that con- 
siderable literature on the dynamics of hologram formation in photorefractive 
materials exists, and that much of it employs more sophisticated modeling than 
that described here. We cite as representative work by Gaylord and colleagues 
at Georgia Tech^^^’^^^ and work by Young and colleagues at the University of 
British Columbia This is in addition to the research cited earlier 
done by Tlttel and colleagues at Rice University . Other authors have also 

noted the exchange of energy between writing beams during hologram formation in 

(4 17 - 19 ) 

photorefractive materials. ’ One group has noted the potential of holo- 

graphic self addition and self-subtraction for use in optical image processing 
applications 
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III. EXPERIMENTAL METHODS FOR THE STUDY OF SELF-SUBTRACTION 


This section concerns the experimental methods used in the program. 
Actual results are described in detail in Sec. IV. We now limit discussion 
to the types of samples studied, the geometry of the experiments, the pro- 
cedures and problems encountered in taking data, and sample preparation 
techniques. 


THE USE OF OPTICAL WAVEGUIDE SAMPLES 


Perhaps the most significant procedural difference between this 

program and other programs that have dealt with photorefractive properties 

of LiNb 03 is the use of optical waveguide samples rather than bulk 

samples. One of several good reasons for using waveguides is that 

potential applications of self-subtraction are attractive to consider imple- 

(3) 

menting in an integrated optical format. Using this format, one has the 
advantages of achieving a high beam intensity using low-optical powers and 
also the advantages of writing holograms that are up to 15-mm thick in sturdy 
yet thin slabs that need be polished on only one side. These features allow 
for the formation of efficient holograms using 0.633-pm recording wavelength, 
which is non-optimum, in writing times of about two minutes. Additionally, 
heat treatments to modify the photorefractive characteristics of the wave- 
guide do not require the tens and even hundreds of hours of thermal equill- 
bration time that thick samples for conventional volume holography do.'' '' 

All these features result from the fact that in the integrated optical for- 
mat only the top few microns of material are exposed to the waveguided beam. 
An assumption that we make is that the nature of photorefractivity is funda- 
mentally the same within the top few microns of material as it is everywhere 
within the sample volume. Although this is an assumption that should be 
checked, we have been comfortable with it in view of the fact that we observe 
similar characteristics for holograms in both integrated optical and 
conventional 3-D optical formats. 
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THE SAMPLE GEOMETRY 


Figure III-l shows the sample-optical beam configuration typical 
of that employed in our experiments. Two coherent beams of 0.633-pm wave- 
length radiation from a He-Ne laser are coupled into the LiNbO^ waveguide 
using a 10-mm wide coupling prism. Interference fringes are established in 
the overlap region which run parallel to the line bisecting the guided beams 
and perpendicular to the optical axis. The light is TE polarized; that is, 
the electric vector is in the plane of the waveguide. Usually the light 
is in the TE^ waveguide mode. It is important that both beams occupy the 
same mode to establish a high-contrast interference pattern and thereby to 
increase the rate of hologram formation. 

The signal beam is the narrower of the two beams shown in Fig. III-l. 
It is typically formed by aperturlng to 1-mm width a beam that is an attenu- 
ated mirror image of the broad reference beam used in the experiment. This 
geometry maximizes the path length and the spatial uniformity of the hologram. 
Moreover, it minimizes the effects associated with the diamond-shaped holo- 
gram region that results from the use of finite beams. Diamond -shaped holo- 
grams exhibit recording characteristics that have special features associated 

(20 21 ) 

with their geometry. * This would complicate the interpretation of our 

results so we minimize the effects by using the wide-beam narrow-beam 
geometry of Fig. III-l. 

In order to use the theory of self-subtraction to analyze our 
experimental data and thereby obtain values for photorefractive parameters, 
it is necessary to know the Intensity of light in each waveguided beam. This 
is a sltxiation in which the integrated optical format presents a problem not 
encountered with the conventional bulk optical format. Generally the power 
throughput of a prism-input and-output coupled waveguide is in the range 
1% - 10%. The power in the waveguide is indicated by the input coupling 
efficiency alone. This could be as small as the throughput itself or as 
large, approximately, as the square root of the throughput. Even if the 
power in the waveguide were known accurately, the intensity varies spatially 
with depth from the waveguide surface according to the mode profile, and an 
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^pp]^opiricit6 civBiTflgG xntsnsity wouXd I13.V6 to bo coiuputod to us6 in tbo f orniuiss 
of the previous section. The net effect is a rather large uncertainty when 
our experimental results are applied to the estimation of parameters in the 
theory. For this reason we have not made parameter determination an important 
part of our experiments. Rather, our interest has been in the shape of the 
curves describing signal— beam output intensity as a function of time, and in a 
comparison of experimental and theoretical curve shapes. The curves scale 

with optical intensity but, in the absence of nonlinear effects, do not change 

( 2 ) 

in shape with changes in optical intensity. This was verified to our 
satisfaction by an experiment conducted during the program. 

For those who may wish to reproduce our experiments, we note that 
typical values of input power in the approximately 4-mm diameter reference 
beam were 5 mW measured before the input coupling prism and 20 pW measured 
after the output coupling prism. Taking the latter number as a lower 
estimate of power in the waveguide and employing an approximate beam cross 

section of 4 mm x 5 pm, the average intensity of the waveguided beam is 

2 

found to be O.lW/cm . This is several orders of magnitude below the inten- 

( 22 ') 

sitles observed to cause significant nonlinear photorefractive effects. ^ ' 

One final aspect of the sample geometry shown in Fig. III-l has to 
do with the period of the hologram fringes, called A, and its relation to the 
Bragg angle, 6 . The relevant formula is 

A = A /2 sin0„ III-l 

O D 

where X is the free-space wavelength of the writing light, and is 

O n 

measured external to the sample. The same formula would apply if 9,, were 

measured in the waveguide sample, except that X^ would then be replaced by 

X, the wavelength in the material. The changes in 9 and X caused by re- 

B O 

fraction exactly cancel each other with the result that an identical value 

for A is calculated. In order to be specific, we take 9_ to be the Bragg 

B 

angle measured external to the waveguide, prior to input coupling. For a 
typical value 9„=15‘’, the hologram period is A = 1.22 pm. 

D 
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THE ARRANGEMENT OF EXPERIMENTAL APPARATUS 


Figure III-2 shows the overall arrangement of the self-subtraction 
experiment, A variable beam splitter was used to divide a collimated laser 
beam of approximately 4 mm diameter into signal and reference beams. The 
beam intensity ratio was in the range from 100:1 to 1:1, with the signal 
beam identified as the weaker of the two. Four mirrors were used to re- 
combine the beams at the LiNbO^ sample as shown in the figure. The additional 
two mirrors over the minimum two required facilitated the adjustments needed 
to couple both beams into the waveguide at once. One mirror was mounted on 
a piezoelectrlcally translated stage to accomplish a voltage controlled 
phase shift. This could then be used to test the sensitivity of a self- 
subtraction hologram to small phase shifts such as those that might some- 
day contain signal information in a practical device. Alternatively, a phase 
shift equal to ir could be applied to peak the Intensity of the subtracted 
beam. The ratio of the peaked signal to the previous subtracted signal is 
called the contrast ratio and is a figure of merit of the subtraction process. 
Moreover, the sudden application of a phase shift is used as the initial 
step in evaluating the ability of the sample to recycle; that is, the ability 
to record new self-subtraction holograms on top of old ones without the need 
for optical or thermal erasure. 

A vertical 1-mm wide slit was placed in the signal beam before 
coupling into the waveguide to reduce the width of the beam without changing 
its intensity. Beam ratios were measured by use of a power meter placed in 
the signal beam prior to aperturing. Consequently, these ratios refer to 
relative intensity in the waveguide and not relative power. The slit was 
found to block approximately 70% of the incident power from the signal beam. 

DATA ACQUISITION 

The typical approach to data acquisition that we employed is de- 
scribed as follows: The output signal beam was directed to a UDT power meter 

for detection. The analogue output from the power meter was used to drive 


34 




FIGURE III-2 


He-Ne LASER (^50 mW) 



ARRANGEMENT FOR SELF-SUBTRACTION STUDIES 





the y-axis of an x-y recorder while the x-axls position was changed at a 
constant rate conmensurate with the hologram- formation time constant. This 
produces the desired permanent record of the evolution of the hologram. The 
method is attractive in its simplicity, yet there are several sources of 
optical errors that have to be taken into account. 

First of all, scattering in the waveguide causes a noticeable 
fraction of the output reference and signal beams to fan out. The strength 
of the scattered component is further enhanced by photorefraction. The net 
effect is that a portion of the signal beam misses the finite collection 
aperture of the power meter and is not detected. An associated problem is 
that a scattered portion of the reference beam can reach the detector and 
be recorded as signal. We have ignored the first problem, as is acceptable 
provided that the fraction of scattered signal power is a constant in time. 
The second problem is more severe because it prevents us from measuring 
unscattered signal power below a certain background level. Since the 
minimum signal power obtained during self-subtraction is a measure of the 
quality of the effect, the background scattering should be eliminated as 
much as possible, and any residual scattering should be measured to be sure 
that it is not a factor in the experiment. 

Generally, background scattering is reduced by moving the detector 
further from the waveguide. Residual levels of background scattering are 
conveniently measured by either of several methods. In one, the signal beam 
is blocked and the scattering from the reference beam is measured as back- 
ground. This method fails if a hologram exists in the sample when the 
measurement is made. In this case light diffracted by the hologram is 
incident on the detector. Since this is not scattered background light, an 
erroneous reading is obtained. In addition the technique suffers from 
the fact that background scattering from the signal beam is not recorded, 
although this component would be negligible in most experiments where the 
beam ratio is on the order of 10:1 or 100:1. 

A second method for measuring background scattering is valid when 
both reference and signal beams illuminate the sample and even when a holo- 
gram has been recorded. In this method, the detector is simply raised and 
lowered relative to its initial position in the path of the signal beam. 
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The average reading in the raised and lowered positions is an indication of 
the background level at the Intermediate position occupied by the signal beam. 
One of the two types of background-level measurements we have discussed 
should be made in all experiments in which the minimum power in the subtracted 
beam is a quantity of interest. 

In addition to the problem of background scattering, another source 
of experimental error to be considered is air turbulence. This effect 
causes fluctuations in the optical lengths of the reference and signal beam 
paths. A legitimate self-subtraction hologram cannot be recorded unless the 
corresponding phase fluctuations are much less than n. To achieve this, the 
experiment of Fig. III-2 was housed beneath a rectangular frame enclosure 
covered with transparent polyethylene sheet. 

A PHASE-TRACKING EXPERIMENT 


Figure III-3 shows a modification of the basic experimental set 
up that was devised to track phase changes of the signal beam as a self- 
subtraction hologram was recorded. This was done as a test of the theory 
of Sec. II. Figures II-6 and 7 of that section indicate that as the sub- 
tracted beam Intensity approaches zero, the temporal rate of change of the 
phase of that beam becomes large. Measurements of the self-subtraction 
intensity alone contain no phase information, so the modification shown in 
Fig. III-3 was established to gather this information. We desired to 
measure any rapidly occuring phase changes that might occur as the self- 
subtracted intensity became small. 

Figure III-3 shows that part of the reference beam is redirected 
by a beam splitter and mixed with a portion of the output signal beam. A 
photomultiplier is set up to detect light from the resulting fringe pattern. 
An aperture is used to limit the field of view to a width small compared to 
the fringe spacing. In this way, signal-beam phase changes are converted to 
amplitude changes detected by the photomultiplier. 

A complication arises from the fact that the signal-beam amplitude 
decreases during subtraction. This will cause the photomultiplier output 
to change even if the phase of the subtracted beam is stable. To avoid this 
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FIGURE III-3. EXPERIMENTAL CONFIGURATION FOR TRACKING 
THE PHASE OF SIGNAL BEAM DURING HOLOGRAM 
RECORDING. 
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complication, a piezoelectrically-translated beam-splitter is used to pro- 
vide the reference beam component to be mixed with the signal beam in this 
experiment. The position of the beam splitter is adjusted to keep a dark 
fringe (or bright fringe) centered on the detector. At any time, the volt- 
age applied to the piezoelectric drive to move the beam splitter the required 
amount is proportional to the accumulated signal -beam phase change. 

A feedback electrical system for tracking the moving fringe pattern 
could no doubt be devised, but for our preliminary investigation, manual 
tracking was acceptable. Details of the experiment results are discussed in 
Sec. IV. 


SAMPLE- PREPARATION TECHNIQUES 

Two types of LiNbO- waveguides were used for the investigations of 

(23) 

this program. These were Ti-indiffused waveguides and Li„0 outdiffused 

(24) ^ 

waveguides . The Ti-dif fused waveguides have grown more popular for 

applications in recent years, owing to the ease with which channel wave- 

(25) 

guides can be formed using standard lithographic techniques , owing to 

the confinement of light closer to the waveguide surface, and owing to the 

existance of TM modes and modes that propagate parallel to the optical axis 

f 26 ^ 

of the LlNbO^ crystal.'' ' These are all features absent in outdiffused 
waveguides, resulting in somewhat restricted use of these waveguides in 
application-oriented programs . While we performed several experiments with 
Ti-dif fused waveguides, the results generally showed oscillatory behavior 
with an overall trend toward holographic self-addition rather than self- 
subtraction. Consequently, we did not pursue this waveguide variety in 
detail. Our feeling was that diffused Ti was participating in the photo- 

refractive process in ways too complicated to effectively analyze using the 

f27 28^ 

results of the limited set of experiments we were performing. ’ ’ We 

emphasized outdiffused samples in the program primarily because they were 
the ones that best exhibited the self-subtraction phenomena of interest. 
Understanding the self-subtraction process in these samples was perceived 
to be the appropriate first step in understanding the process in waveguides 
generally. 
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Those wishing to repeat our experiments may note the waveguide-fabrication 
methods employed. Ti diffused waveguides are formed by e-beam evaporation 

O 

of about 200 A Ti on the LiNb 03 substrate surface and subsequent diffusion 
of the metal into the substrate at 950° or 1000° C for about 3 h. During 
the diffusion, O 2 gas is passed over the sample. After diffusion the 
sample is quenched to about 600° C in an interval of approximately 15 sec. 

It is then allowed to cool to room temperature at its natural rate or at an 
accelerated rate requiring as little as 30 min. Usually O 2 gas remains 
flowing over the sample during the cool down. Although the flowing gas has 
no effect on waveguide characteristics at temperatures of 600° C or less, it 
can have significant effect on the photorefractive characteristics of a 
sample. We have demonstrated this by subjecting waveguides to reducing 
treatments in Ar and oxidizing treatments in flowing O 2 at temperatures up 
to 600° C and for times up to 5 h. 

The motive for conducting these heat treatments is discussed in 

3+ 

Ref. 8. A reducing heat treatment is believed to add electrons to Fe iron 

2 + 

traps in the LiNbO^, thus converting them to Fe donors. This increases 
the photorefractive sensitivity of the crystal at the same time it decreases 
the supply of available Fe electron traps. In the simplest model, the 
mean free path L of the photo-excited electron varies Inversely with the 
number of available traps . We thus have a mechanism by which a reducing 
thermal treatment can increase L. According to the theory of Sec. II, 

L/A>>1 is a condition desired for the best self-subtraction. We thus expect 
a reducing heat treatment in flowing Ar to Improve a sample’s self-subtraction 
characteristics. A counter line of reasoning is used to show that an 
oxidizing thermal treatment can reduce L and hence degrade a sample's 
capacity for self-subtraction. 

Although we employed a number of oxidizing and reducing thermal 
treatments on several LiNbO^ waveguides, in only one case were we successful 
in improving self-subtraction performance. In general our results were not 
consistent enough to permit the identification of trends to use for guidance 
in specifying sample preparation methods. The details of particular ex- 
periments, however, are provided in the following section. 
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IV. RESULTS OF EXPERIMENTAL STUDIES OF SELF-SUBTRACTION 


This section contains the principal experimental results of the 
program. Approximately 50% of the work was done on a particular outdiffused 
waveguide sample, #96, which exhibited the best self-subtraction of any 
sample studied. We begin the section by describing our detailed results 
with this waveguide. Including its history, its potential for self-subtraction 
data processing, the effects of various experimental configurations on per- 
formance, and effects of various heat treatments on performance. Then we 
will compare our results with sample #96 to results obtained using similarly 
fabricated samples to show the extent of reproducibility and predictability 
of self-subtraction phenomena. We will briefly describe the nature of our 
results with a Ti-lndiffused LiNbO^ waveguide, and we will close by relating 
our experimental results to the theoretical predictions of Sec. II. 

EXPERIMENTAL STUDIES OF SELF- SUBTRACTION IN 
OUTDIFFUSED WAVEGUIDE SAMPLE #96 

Sample #96, purchased "off the shelf" from Crystal Technology, is 
a 25.4 mmx 12.7 mmx 1mm y-cut slab that was probably intented for acoustic 
applications. The polish was Inferior to the optical quality polish we 
receive on specially-prescribed crystals for integrated optical applications. 
In spite of those unpromising origins, #96 exhibited the best self-subtraction 
of any sample we studied. The waveguide was made by thermal outdlf fusion at 
950° C for 1 h in flowing O 2 . The sample was quenched to 600° C and then 
allowed to cool slowly to room temperature, probably in flowing O 2 at all 
times . 

The first hologram recorded in the sample revealed a slow, 
oscillatory signal-intensity-vs-time curve that tended slowly toward the 
null associated with self-subtraction. A second run with a reference-to- 
signal beam ratio of 7.5/1 (down from 33/1) showed a barely noticeable 
oscillation, and the time required for the signal-beam intensity to fall to 
1/e of its initial value was about 320 sec. Following these runs, the sample 
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was subjected to a 2h anneal in flowing Ar at 500° C to see of any improve- 

/Q\ 

ment in self-subtraction performance could be obtained . ^ ^ In fact, the 
oscillatory characteristics returned stronger than previously. However, 
after six days of sample inactivity, a series of new runs were made that 
showed little or no oscillation and good quality self-subtraction. These 
runs are the subject of the following subsection. 

Tests of Sample #96 as a Self-Subtraction Processor 

We made a series of measurements with sample #96 in which the 
crystal performed the operations expected of a self-subtraction processor. 

That is, repeated self-subtraction holograms were formed after the application 
of repeated phase shifts and optical erasures. Performance was characterized 
by measurements of (1) the time constant for subtraction, (2) the contrast 
ratio, indicating the extent of subtraction attained, (3) the stability of 
a self- subtraction hologram over long periods of constant writing and read- 
out, and (4) the reproducibility of the subtraction phenomena after multiple 
optical erasures and phase shifts. 

Most of the results to be described were obtained during the course 
of one week of investigation. The activities of that week are summarized 
by Table IV-1. In the context of the table, "erase" means optical erasure 
with the signal beam blocked. "Rewrite" means writing a new self-subtraction 
hologram on top of an old one following the application of a phase shift to 
one of the writing beams. "Erase? (Tesla)" means that a Tesla coll was 
used to attempt erasure of a waveguide hologram. We did not attempt to 
verify that the hologram was erased, but we did notice an increase in re- 
cording sensitivity after the experiment that could possibly have resulted 
from the electron bombardment of the surface. 

w 

Time Constants, Contrast Ratios, and Reproducibility of Self-Subtraction 
Holograms 


Figure IV-1 shows the formation of the self-subtraction hologram 
associated with the second rewrite of day 7, immediately following an applied 
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TABLE IV-1. SELF- SUBTRACTION EXPERIMENTAL ACTIVITY USING LiNbO„ 
OUTDIFFUSED WAVEGUIDE SAMPLE //96 ^ 


Day 

1 

2 

3 

4 and 5 

6 


7 

8 

Activity 

Write 

Erase 

Write 

Inactive 

Rewrite 

Rewrite 

(? - shift) 

Erase 


Erase 

Write 

Erase 


Erase 

Rewrite 

(it - shift) 

Write 


Write 

Erase 

Write 


Write 

Rewrite 

(tt/4 - shift) 





Erase? (Tesla) 



Rewrite 

(n - shift) 





VJrite 



Rewrite 

(= Beam Intensities) 



U) 



SIGNAL-BEAM POWER, yW 



FIGURE IV-1 . POWER-VS-TIME FOR A SUBTRACTED BEAM THAT EXPERIENCES A ir-PHASE SHIFT • 
AT t = 0 AND WRITES TO A RENEWED CONDITION OF SELF SUBTRACTION. 


phase shift of ir to peak the signal-beam intensity. The duration of writing 
encompassed four successive sweeps of the recorder, with the last sweep 
interupted for about 30 min. Just prior to application of the phase shift, 
the previous hologram had achieved a subtraction level of about 0.06 yW. 
Immediately after the phase shift the signal beam intensity was 2.6 yW, 
indicating a contrast ratio for the initial hologram of 2.6/0.06 = 43.7. 

Based on subsequent experiments, it appears possible that there was back- 
ground scattering contribution of 0.01 yW to the minimum measured signal of 
0.06 yW. In this case the contrast ratio would be as large as 2.6/0.05 = 52.0. 

The signal for the rewritten hologram of Fig. IV-1 falls to 1/e of its 
peak value after 144 sec. The signal is 0.1 yW after 946 sec., and after 
50 min. of recording it is reduced to 0.05 yW. As mentioned above, as much 
as 0.01 yW of this subtracted signal could have been associated with back- 
ground noise. No tendency toward instability is noted after 50 min. of 
continuous recording. 

Following this experiment a phase shift of about ir/4 was applied to 
one of the writing beams, and a third self-subtraction hologram was formed. 

The writing characteristics for this hologram were similar to those of 
Fig. IV-1, despite the smaller value of the applied phase shift. The sub- 
tracted power was a maximum of 0.06 yW, and no tendency toward instability 
was noted after more than 30 min recording time. The time constant for 
recording this hologram was 174 sec., and the contrast ratio at the end of 
recording was a minimum of 40.0. 

Finally another tt phase shift was applied and yet a fourth self- 
subtraction hologram was recorded. The time constant was 139 sec. and the 
final subtracted power was 0.06 yW, as before. The signal power-versus- 
time curves for holograms three and four of day seven are shown in Fig. IV-2 
and 3, respectively. 

These results were obtained with writing beams having unequal 

intensities in the approximate ratio 14:1. The more intense reference beam 

was 4 mm wide and the signal beam was 1 mm wide. The power in the output 

reference beam was 20.8 yW, from which we estimate the waveguide intensity 

2 

to be on the order of O.lW/cm . The external Bragg angle of the hologram 
was 14°, corresponding to a fringe spacing of 1.3 ym. The beam geometry is 
that shoxjn in Fig. III-l. 
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IHAT EXPERIENCES A ir/4 


RENEWED CONDITION OF 



signal-beam power, yW 



FIGURE IV-3. POWER-VS-TIME FOR A SUBTRACTED BEAM THAT EXPERIENCES A 
SECOND TT PHASE SHIFT AND WRITES TO A RENEWED CONDITION 
OF SELF SUBTRACTION. 
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The last hologram recorded on day 7 and the only hologram of day 8 
differed from the above holograms in that nearly equal writing beam intensities 
were employed. A slight tendency toward instability was noted. To be 
specific, in the final run a contrast ratio close to 57 was reduced to 27 
by continued writing at which level stabilization appeared to have been 
achieved . 

Earlier experiments, namely, those occurring on days 1-5 of 
Table IV-1, differed from later experiments in that equal widths for signal 
and reference beams were employed. The widths, were probably close to 3 ram. 

The overall area of the hologram is several times larger as a result of the 
use of two wide beams, and this appears to be associated with a slight 
reduction in the self-subtraction time constant. With reference and signal 
beam intensities equal to each other but less than the intensities used in 
the last runs described above, the time constant was reduced from 140 sec. 
to about 100 sec. Further reduction occured with the use of unequal beam 
intensities. The shortest time constant of about 60 sec was obtained with 
a 73/1 beam ratio following Tesla-coil exposure of the waveguide surface. 

We cannot speculate as to whether the large beam ratio or the exposure 
dominated in producing the reduced time constant, but this is a question 
worthy of future investigation. Time constants lower than tens of seconds 
will probably be desirable for most applications. 

Sensitivity of Self-Sub traction Holograms to Applied Phase Shifts 

Figure IV-4 shows the results of an experiment to test the 
sensitivity of a self-subtraction hologram to applied phase shifts. The 
hologram is that associated with Fig. IV-3. The initial subtracted intensity 
is about 0.06 pW. At intervals averaging 40 sec. an external phase shift 
is applied to one of the writing beams and then removed. Each phase shift 
is larger than the previous one by a certain increment A<() which can be 
determined from the data. We assume that the intensity of the n-th pulse 
in Fig. IV-4 is described by 

I = + 0.51 ^ [1 - cos (6 + nA<^)], (IV-1) 

n bg max o 
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01234 56 7 

TIME, min 

FIGURE IV-4. INTENSITY VS TIME FOR A SUBTRACTED BEAM IN THE PRESENCE 
OF EXTERNALLY APPLIED PHASE SHIFTS. 


VO 



where is the background Intensity prior to the phase shift, is the 

peak signal associated with a it phase shift, is any phase affset associ- 
ated with the self-subtraction mechanism, and Acfi is the phase increment. 

Figure IV-5 shows a plot of (I - L )/l versus the phase index 
n. The theoretical fit is obtained using = 0 and Atj) = 22° in Eq. (1) . 
There is fairly good agreement at lower values of nA(}>. From Fig. IV-4, note 
that the first phase shift A<j> = 22° results in a signal comparable to the 
background intensity level. We have, therefore S/N = 1 when A(j) = 22°. 

This describes the sensitivity of the sample as a phase-shift processor. 

Figure IV-5 shows increasing discrepancy between theory and ex- 
periment at values of nAi}> greater than 180°. This may indicate a slight 
degradation of the hologram caused by the repeated phase shifts. It is 
interesting to note that following the last phase shift the minimum sub- 
tracted signal is observed at =9° rather than <|) = 0°, as is the case 

o o 

initially. That is, an offset in the relative phase required for self- 
subtraction develops upon repeated applied phase shifts. If one allows the 
beams to continue writing, however, the intensity at <J>q = 9° increases at 
first and then falls back to its minimum value where it remains stable. 


Stability of a Self-Subtraction Hologram to Continuous Seven-Hour Readout. 

The experiments summarized by Fig. IV-1,2, and 3 show that self 
subtraction holograms are stable under continuous readout for times up to 
65 min. In applications, self-subtraction processors could very possibly 
be subjected to continuous readout over a period of months, during which 
they should remain stable. By use of the word stable we mean to infer that 
no significant deviations from the minimum subtracted signal occur except 
those which are triggered by phase shifts in the writing beams. 

We investigated the stability of a self-subtraction hologram in 
sample #96 for a maximum interval of seven hours. Figure IV-6 shows the 
first 30 min. of hologram recording, during which time a subtracted intensity 
level of 0.04 pW is achieved. Figure IV-7 then shows the subtracted in- 
tensity-versus-tlme variation for the next 6 h 30 min. The average sub- 
tracted signal of 0.03 pW is stable to within + 15%. The degree of 
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POWER-VS-TIME CHARACTERISTICS OBSERVED DURING THE 
FIRST 30 MIN OF A 7 H STABILITY TEST. 





TIME , hours 



2. 5-3. 5 




3. 5-4. 5 



5. 5-6. 5 

6. 5-7.0 
T7 PHASE SHIFT 



FIGURE IV-7. SUBTRACTED POWER (VERTICAL) VS TIME (HORIZONTAL) IN 
THE INTERVAL 0. 5-7.0 h. AVERAGE SUBTRACTED POWER 
= 0.03 yW. 
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fluctuation about the average signal changes in time and is apparently the 
result of activity in the laboratory, such as people working, or changes in 
the environment associated with the air-conditioning system. At the end of . 
the run a tt/ 2 phase shift was applied and a contrast ratio of 43 measured. 

This is comparable to the contrast ratio achieved after much shorter runs. 

The beam ratio used to make this run was 14/1, as measured at the waveguide 
output. A value 34/1 was measured at the waveguide input. The difference 
indicates nonuniform coupling of the reference and signal beams. This is not 
Important to the experiment but indicates a large uncertainty as to the cor- 
rect beam ratio in the waveguide. Other parameters of the experiment were 
the same as_ those for the experiments summarized by Figs. IV-1,2, and 3. 

Sensitivity of Self-Subtraction Holography to Bragg Angle 

The results described to this point suggest that self-subtraction 
characteristics are rather insensitive to parameters such as beam width and 
beam ratio. We now begin a discussion of parameters to which self-subtraction 
is fairly sensitive, beginning with the Bragg angle used in recording; that 
is, the half angle between the two intersecting beams as measured in air. 

The motive for studying the dependence of self-subtraction on the 
Bragg angle, 6 , is the important role played by the quantity tan ^ KL in 

D 

the theory of Sec. II, where 

KL = (4t:L/A^) sinB^, (IV-2) 

L is the mean free path of conduction electrons and is the free-space 

optical wavelength. Values of KL larger than one are most compatible with 

good self-subtraction, at least according to the photorefractive model we 

employed. One can vary KL by changing the mean path length L or the magnitude 

of the grating wavevector, K = (4r/X ) sin0 . An easy way to vary the latter 

parameter is to change the Bragg angle, 0„. Since our initial choice 0„ = 14“ 

B B 

gave good self subtraction, indicating KL3^1, we thought that 10 X reductions 

in 0g might be necessary to reduce tan ^ KL enough to see an effect. However, 

we found that significant variations in self-subtraction characteristics 

occured with only 20% variations in 0 . 

B 
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We tested five different Bragg angles in the range 5.5° to 17.9°. 
For = 17.9° and 13.9°, the subtracted beam intensity showed a generally 
monotonic decrease toward the saturated level. For 11.1°, 7.25° and 5.5°, 
a short-lived initial decrease in the subtracted intensity was followed by 
a large increase (holographic self-addition) and then the subtraction 
process took over. The variation in intensity as a function of writing time 
for these experiments is shown in Fig. IV-8. The photoref ractive model of 
Sec. II does not appear to provide a basis for interpreting these results. 

Sensitivity of Self-Subtraction Holography to Thermal Oxidizing and 
Reducing Treatments. 

The variation of 0 was used, as described above, to accomplish a 
B 

variation in the photoref ractive parameter KL of Eq. 2 by causing a change 
in K. Oxidizing and reducing heat treatments can change the product KL 
through their effect on the electron path length L. We therefore were 
motivated to study the variations in self-subtraction caused by oxidizing 
and reducing thermal treatment of sample //96 and other samples in oxygen 
and argon, respectively. The details of these thermal treatments were 
described in Sec. III. 

In this section we have already indicated that a reducing treatment 
performed on waveguide #96 after it was fabricated helped to improve its 
self-subtraction characteristics. We now describe the results of other 
thermal treatments on #96 that took place in the time frame 4.5 - 5.0 months 
after the last experiments described above, those summarized by Fig. IV-8. 

The photorefractive properties of the sample do not appear to have degraded 
over the 4.5 month delay between experiments. However, a 5 h heat treatment 
in flowing Ar at 600° C was found to cause a serious impairment in the self- 
subtraction characteristics of the sample. With a Bragg angle of 18°, the 
initial reduction of signal beam intensity with recording time progressed at 
a rate slower than previously and reversed itself before 50% of the initial 
signal level was reached. No subsequent tendency toward subtraction was 
observed during the observation time, which was approximately twice the time 
that the initial signal intensity decreased. In several runs the final 
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FIGURE IV-8. SELF-SUBTRACTION CHARACTERISTICS OF HOLOGRAMS IN OUTDIFFUSED LiNb03 WAVEGUIDE 
§ 96 , MEASURED FOR FIVE DIFFERENT BRAGG ANGLES. 


intensicy observed was greater than the initial intensity (holographic self 
addition) . It is possible that this situation would have reversed if we had 
taken data for a long enough period. Figure IV-8 provides evidence for this' 
possibility. However, the time required would have been so great as to 
render the sample useless as a self-subtraction processor, so we did not 
continue the measurement. 

When a second heat treatment of 5 h in flowing Ar at 600° C was 

given to sample #96, the self- sub traction characteristics of the sample did 

not appear to change significantly from the degraded characteristics 

observed following the first reducing heat treatment. As described in 

Sec. Ill, it was anticipated that these reducing heat treatments would have 

3+ 2+ 

a beneficial effect on self-subtraction by converting Fe traps to Fe 
donors, thus increasing both photorefractive sensitivity and path length. 
Apparently this simple-minded viewpoint is not sufficient to describe the 
actual dynamics of the photorefractive process. Subsequently, we provided 
sample #96 with an oxidizing heat treatment to see if we could reverse the 
self-subtraction degradation caused by the two reducing heat treatments. 

The oxidation was carried out over a period of 4 h, with the sample main- 
tained in flowing O 2 at 600° C. The effect of this treatment was, if any- 
thing, to reduce the exchange of energy between writing beams. No tendency 
toward either self subtraction or self addition was observed. 

Our final experiment with sample #96 was to reproduce the conditions 
used in making the initial waveguide; that is, the sample was heated in 
flowing O 2 at 950° C. The length of the heat treatment was 1/2 h, in contrast 
to the 1 h treatment used to make the initial waveguide. The sample was pre- 
heated in O 2 at 600° C for 2.75 h prior to raising the temperature to 950° C. 
At the end of the heat treatment the sample was rapidly cooled to 600° C and 
then brought to room temperature over a period of 1/2 h. Oxygen gas remained 
flowing over the sample throughout the process. Only one self-sub traction 
run for the sample was performed following this treatment and that run 
indicated little or no tendency toward self-subtraction. 

In summary, a variety of thermal treatments were performed on 
sample #96. A 2 h reduction in flowing Ar at 500° C was performed after 
the outdiffused waveguide was fabricated, and this was the only beneficial 
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heat treatment from the point of view of optimizing self-subtraction. Even 
in this case, beneficial effects appeared only after a 6 day time delay, and 
it is possible that the thermal treatment had little to do with their 
appearance. 

Although the good self-subtraction characteristics of sample #96 
remained in force for 4.5 months, the first subsequent thermal anneal at 
600® C removed them, and no thermal anneal attempted thereafter was able 
to restore them. In fact, no significant trends were observed following 
successive oxidizing or reducing treatments. As a result we can not 
confidently predict the most suitable fabrication method for a waveguide 
to be used for self-subtraction, nor can we state the optimum self-subtraction 
performance to be expected from LlNbO^ waveguides. However, we did work 
with other waveguide samples than #96, and these tended to show fair self- 
subtraction in their as-fabricated state when they were fabricated by methods 
similar to those used for sample #96. This provides us with a starting 
point for any future research that may be undertaken on self-subtraction. 

It is suggested that thermal heat treatments of shorter duration and lower 
temperatures be employed to insure that dramatic changes in self-sub traction 
performance are avoided. In this way, trends may be identified and an 
optimized set of fabrication conditions may result. In the following sub- 
section we present, for guidance, the results of our experiments on samples 
other than #96. 


EXPERIMENTAL STUDIES OF SELF-SUBTRACTION IN 
OUTDIFFUSED WAVEGUIDE SAMPLES #97, 100, AND 162 


Performance of Sample #97 


Sample #97 was outdlffused using the same procedure employed for 
sample #96; namely, a 950° C anneal in flowing O 2 for 1 h, followed by a 
rapid quench to 600® C and a slow cool to room temperature. The light anneal 
in Ar that was applied to #96 after its fabrication was not applied to #97. 
Perhaps as a result, the performance of #97 never matched the best perfor- 
mance of #96. Instead, #97 behaved qualitatively similar to the way #96 


58 



behaved in the as-fabricated state, showing slow subtraction with a tenden- 
cy toward oscillation or instability. Self-subtraction characteristics ap- 
parently did not benefit from a 6.6 month delay between fabrication and 
testing. 

The ratio of initial signal-beam intensity to reference beam in- 
tensity that gave the best subtraction was 1/10. In this case the self- 
subtraction time constant was about 10 min. However, there appeared to be 
a tendency toward instability since the signal-beam intensity was climbing 
slowly when the run was terminated at 27 min. In the case of 1/1 or 1/100 
beam ratios, the steady-state trend seemed to be toward a constant, unsub- 
tracted signal beam. We recall that sample #96 did not show a similar 
dependence of subtraction capability on the beam ratio. 

Following these intial measurements, sample #97 was subjected to 
two 5 h heat treatments in Ar at 600° C. These did not result in any 
significant change in the self-subtraction characteristics of the sample. 
However, a third heat treatment, in flowing oxygen at 600° C for 4 h, re- 
sulted in the very dramatic changes shown in Fig. IV-9. The signal beam 
intensity-vs-time characteristic is quasi-per iodic. Deviations from exact 
periodicity are fairly regular. Unfortunately, the physical changes in the 
crystal that served to produce this interesting effect remain hidden from 
view. It would, in the future, be useful to perform a detailed analysis of 
the theory of Sec. II to see if effects like those observed here can be 
explained using the simple photorefractive model on which that theory is 
based. Ultimately, it will be desirable to experimentally characterize the 
oxidation-reduction state of the crystalline Impurities to obtain more 
information about the source of behavior like that shown in Fig. IV-9. 

Performance of Sample #100 

Sample #100 was an outdiffused LiNbO^ waveguide formed in "off the 
shelf" material similar to samples #96 and 97. The outdif fusion heat treat- 
ment, however, was carried out very differently. The sample was heated to 
1100° C, at which point the oven was turned off. The time spent in going 
from 950° C to 1100° C was 7 min., and the time spent in cooling to 950° C 
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FIGURE IV-9. SIGNAL-BEAM INTENSITY VS TIME FOR SAMPLE #97 
FOLLOWING AN OXIDIZING THERMAL TREATMENT. 




again was 6 min. The waveguide expected to result from this procedure is 

(29) 

shallower than the conventional outdif fused waveguide. Consequently, 

its photorefractive properties should be more sensitive to the oxidation- 
reduction state of impurities at the waveguide surface. 

Several runs were made using this sample. One run was made using 
the waveguide as fabricated; another was made after a four hour oxidation 
carried out at 600° C. The most notable feature of the signal-intensity- 
vs-time characteristics of holograms formed in this waveguide was their 
lack of consistency. In some runs the signal intensity decreased in time 
as expected for self-subtraction. When it did so, however, two time con- 
stants were in evidence for different rims; one on the order of minutes, 
the other on the order of hours. In other hologram-formation experiments, 
the signal-intensity versus time curve showed oscillatory behavior with a 
characteristic period on the order of 20 min. The oscillations were not, 
however, periodic or quasi-periodlc like those of Fig. IV-9. 

Because of these inconsistent and unexplicable results, we ceased 
experimenting with sample #100. One conclusion that may be drawn from the 
work is that the waveguide fabrication treatment employed for the sample 
was not conducive to good self -sub traction. 

Performance of Sample #162 


Sample #162 differed from previously discussed samples in that it 

3 

was special-ordered for integrated-optics applications. The 25 x 25 x 3 mm 

3 

substrate was sawed into two similar 12 x 25 x 3 mm substrates, designated 
162A and 162B. The former was subjected to a 1 h outdif fusion treatment 
similar to that used for the preparation of #96 and #97. The latter was 
subjected to a 30 min. outdif fusion to see what effects, if any, might be 
associated with the time of the thermal treatment, other factors remaining 
constant. 

Most experimental data were taken using sample #162A over a period 
of two weeks. Although the range of experiments performed was not as ex- 
haustive as with sample #96, sample 162A appeared to come close to matching 
the performance level of #96. For example, self-subtraction holograms were 


61 



repeatedly written following the application of it phase shifts. The time 
constants for reaching the subtracted condition were typically 2 min. 

Contrast ratios were not measured but appeared to be comparable to those 
obtained using sample #96. 

What is most interesting about sample #162A is that its perfor- 
mance as a self-subtraction waveguide improved after the first three days 
of experimentation with it. During these days the sample exhibited self- 
subtraction, but an occasional tendency of the signal beam to oscillate in 
intensity rather than to subtract was observed. This tendency vanished in 
experiments performed about one week following initial experiments with the 
sample. A similar delay in the onset of good self-subtraction performance 
was also observed in sample #96, though in this case an intervening thermal 
reduction treatment could have had some effect. The possibility is that 
optical preconditioning of the sample is necessary to produce good self- 
subtraction characteristics. Other samples tested in the program were not 
evaluated over the course of a week, at least without complicating inter- 
vening thermal treatments. Consequently, any beneficial effects of optical 
preconditioning were not observed. 

Sample #162B, similar to #162A except for a shorter waveguide 
fabrication time, was tested only over a two day span. It was not foimd to 
exhibit good self-subtraction, even when compared to the initial results for 
#162A. The lesson seems to be that shorter outdlf fusion times are not helpful 
for forming waveguides with good self-subtraction characteristics. 

Fig. IV-10 shows the results of experiments performed using sample 
#162A in which a sequence of five abrupt phase shifts is applied to one of 
the recording beams. Following each phase shift, the signal-beam intensity 
(lower curve in the figure) begins to decrease as required for self-subtraction. 
Also shown in the figure is the reference beam intensity (upper curve) which 
increases as the signal-beam intensity decreases (holographic self addition). 
This is in accord with the requirement that optical energy be conserved be- 
tween the two writing beams. In Fig. IV-10 the vertical scales for the 
reference and signal-beams are not Identical. Consequently the different 
amplitudes for the two beams are not in violation of the energy-conservation 
requirement . 
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TIME, MIN 


FIGURE IV-10. TEMPORAL VARIATION OF SIGNAL AND REFERENCE-BEAM 
INTENSITY IN THE PRESENCE OF FIVE APPLIED PHASE 
SHIFTS. 



EXPERIMENTAL STUDIES OF SELF-SUBTRACTION IN 
Ti-INDIFFUSED WAVEGUIDE SAMPLE #135 


Waveguides formed in LiNbO. by the thermal diffusion of oxidized 

^ (23 25—28) 

Ti films are commonly used in integrated optic applications. ’ ' Two 

good reasons for this are (1) the relative ease with which channel waveguide 
structures are formed using photolithography and (2) the confinement of the 
light close to the waveguide surface where it can be more readily switched 
or modulated by acoustic and electric perturbations. For holographic 
applications, however, the Ti dopant adds a degree of complexity because 

it acts as an impurity in altering the photorefractive properties of the 

(27,28) 

waveguxde . 

The TirLiNbO^ waveguide that we studied in most detail was formed 
by the thermal diffusion of a 270A Ti film at 950° C for 3 h flowing O 2 . We 
looked for self-subtraction phenomena at Bragg angles of 5.5° and 20.6°. We 
took data for both positive and negative c-axis orientations to insure that 
the proper conditions for self-subtraction were achieved at least part of the 
time. Nevertheless, at a 5.5° Bragg angle, holographic self-addition was 
encountered for both c-axis orientations. In each case, the weaker of two 
beams used to write the grating was amplified until a saturation intensity 
was reached. The saturation intensity was 32 times the initial intensity in 
one case, and 9.6 times the initial intensity in the other. 

When a Bragg angle of 20.6° was employed, we observed that the 
intensity of the weaker of the two writing beams oscillated in time. For one 
c-axis orientation, successive oscillations appeared to be damped. For the 
other orientation successive oscillations were as shown in Fig. IV-11. The 
quasi-periodic nature of the oscillations, with increasing amplitude and 
period, is reminiscent of the results obtained with outdiffused waveguide 
#97, shown in Fig. IV-9. 
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COMPARISON OF THEORY AND EXPERIMENT 


Subtracted-Beam Intensity Versus Time 

We have previously indicated that the theory for self-subtraction 
presented in Sec. II has not been exercised enough for us to learn all of 
its ramifications. We do not know, for example, if the theory can describe 
the quasi-periodic oscillations of Figs. IV-9 and IV-11. In so far as the 
theory has been exercised, it predicts that self-subtraction does exist for 
a certain reasonable photorefractive model. This is consistent with the 
fact that we experimentally observe self-subtraction, except that our 
observations show an exponential-like decay to the subtracted state whereas 
the theory predicts a damped-oscillatory decay. Apparently, some feature 
has been left out of the physical model which might account for the observed 
smoothing of the subtraction curve. Only in one experiment did we observe a 
damped oscillatory approach to the subtracted state. This result, shown in 
Fig. IV-12, was obtained using sample #162A. It was not reproduced in sub- 
sequent experiments with that sample. 

Subtracted— Beam Phase Versus Time 


A measurement of subtracted Intensity as a function of time, such 
as that shown in Fig. IV-12, conveys no information about the phase of the 
subtracted beam. A measurement of the phase as a function of time there- 
fore provides a second means of comparison of theory and experiment, one 
which we explored at the close of the program. 

The theoretical results of Fig. II-6 and II-7 shows that the phase 
of the subtracted beam is relatively constant except at times when the sub- 
tracted intensity is near a minimum. Our thought was that a sequence of 
rapid phase variations might occur in our experiments even though, for one 
reason or another, local minima in the subtracted beam intensity were washed 
out. 


66 



RELATIVE SIGNAL 





The measurement of the optical phase of the beam imdergolng sub- 
traction was performed using the experimental set-up shown in Fig. III-3. 

A piezoelectrically positioned mirror was employed to keep a dark inter- 
ference fringe centered on a detector as the fringe pattern attempted to 
move in response to changes in the subtracted-beam phase. The measured 
phase was therefore directly proportional to the voltage applied to the 
piezoelectric positioner. Calibration is performed by observing the 
voltage required to cause the detector to see a bright fringe (it phase shift). 

The result of the experiment is shown in Fig. IV-13. Figure IV-13a 
shows the subtracted beam phase as a function of time on a linear scale 
calibrated in radians, while Fig. IV-13b shows the subtracted-beam intensity 
as a function of time on a relative linear scale. A total phase shift of 
approximately 2ir is observed over the recording interval. Theory predicts 
that a phase shift of this magnitude should result from two segregated it 
phase shifts that occur in the neighborhood of minima in the subtracted- 
beam intensity. In fact a rather continuous change in phase is observed. 

This is not surprising in view of the fact that the subtracted-beam intensity 
variation shows no apparent minima. We may speculate that the same physical 
mechanism that removes the oscillations in intensity also acts to smooth out 
sudden variations in phase. It is, however, somewhat reassuring to find 
that total phase changes as large as 2ir are observed at the same time that 
a significant measure of self-subtraction is observed, since this feature, 
at least, is compatible with the theoretical predictions of Sec. II. 
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V. SUMMARY AND CONCLUSIONS 


This report describes the results of an experimental and theoretical 
program to characterize and understand the holographic process termed self- 
subtraction. The program was very successful in its characterization objec- 
tives. A variety of LiNbO^ waveguide samples were evaluated with respect to 
self-subtraction parameters including (1) time constant for the development of 
a self-subtraction hologram, (2) minimum attainable null of the subtracted- 
beam intensity, as measured by the contrast ratio, (3) number of different 
self-subtraction holograms that can be recorded prior to the onset of material 
fatigue, and (4) stability of a self-subtraction hologram to long-term readout. 

Sample #96 was characterized in depth regarding all four of these 
parameters. Other samples were tested until it was verified that they fell 
short of the performance of #96 in one or more areas. Sample #162A appeared 
comparable to #96, suggesting at least a limited degree of reproducibility, 
but it was not tested in depth. 

The performance of sample #96 was very encouraging with regard to the 

possibility of someday using self-subtraction as a method for discriminating 

(2 3) 

against slowly varying background information in an optical processor ’ . We 

found self-subtraction time constants in the vicinity of 150 sec, contrast 
ratios in excess of 40:1, stability of a self-subtraction hologram during read- 
out for at least 7 h (the longest time examined) , and reproducibility of the 
subtraction phenomena after five optical erasures in a single day (the largest 
number attempted) . While these figures do not guarantee the ultimate utility 
of self subtraction as a "smart" optical-processing technique, we do regard 
their appearance in an initial program to be extremely promising. 

On the negative side of the ledger, the present program had the 
dual objective of understanding the physical basis for self subtraction as well 
as characterizing the effect. It was hoped that this understanding would pro- 
gress to the point where optimization of self-subtraction characteristics would 
be possible. With regard to achieving this objective the results of the program 
were, unfortunately, less than desired. 
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Theoretical calculations based on a simple but realistic photorefrac- 
tive model showed that the beam undergoing self subtraction decreases in inten- 
sity in an oscillatory fashion. Experimentally, the decrease was more nearly 
exponential. Attempts to vary the photorefractive material properties by oxi- 
dizing and reducing heat treatments, and by the use of modified waveguide fab- 
rication techniques, produced major changes in self-subtraction characteristics 
that could not be correlated with the theory. Even simple changes in the 
experimental geometry, such as changing the Bragg angle, produced unanticipated 
major changes in self subtraction. 

Accordingly, we can not state that self subtraction is well understood 
as a result of the present work. I'fhile the effect has been demonstrated to a 
high level of satisfaction, our opinion is that improved understanding of the 
mechanism remains an important factor in any decision one might make as to the 
potential of self subtraction for optical processing. This understanding should 
come in part from experiments like those we have performed, with smaller pertur- 
bations of the crystals between measurements so that trends can be identified. 

In addition some type of measurement of the impurity state of the sample should 
be made with each run so that the effects of various heat treatment are known 
rather than inferred. 

Theoretically, a more sophisticated photorefractive model than that 
employed here should be introduced. Attempts should be made to understand the 
quasiperiodic oscillations that sometimes occur in the signal-beam intensity 
in addition to the self-subtraction effect. Also, attempts to understand the 
physical basis for the observed sensitivity of self subtraction to the Bragg 
angle should be made. 

With this type of effort, one might reasonably expect to optimize the 
self-subtraction effect for use in applications. As an additional benefit, con- 
siderable information will be added to the current understanding of photorefrac- 
tivity, a holographic process that has relevance to many optical storage and 
processing systems, including some that function without the need or benefit of 
self subtraction. 
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